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Abstract: Minimally invasive, specific measurement of cellular energy 
metabolism is crucial for understanding cerebral pathophysiology. Here, we 
present high-resolution, in vivo observations of autofluorescence lifetime as 
a biomarker of cerebral energy metabolism in exposed rat cortices. We 
describe a customized two-photon imaging system with time correlated 
single photon counting detection and specialized software for modeling 
multiple-component fits of fluorescence decay and monitoring their 
transient behaviors. In vivo cerebral NADH fluorescence suggests the 
presence of four distinct components, which respond differently to brief 
periods of anoxia and likely indicate different enzymatic formulations. 
Individual components show potential as indicators of specific molecular 
pathways involved in oxidative metabolism. 
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1. Introduction 

Identifying and characterizing the complex cellular and molecular mechanisms of cerebral 
energy metabolism is critical to better understand brain function and its pathological 
alterations [1]. In order to accurately interpret the noninvasive, macroscopic observations of 
cerebral oxygenation and blood flow obtained by functional neuroimaging techniques such as 
positron emission tomography (PET) and functional magnetic resonance imaging (fMRI), 
more comprehensive insight into metabolic activity at a cellular level is crucial. Additionally, 
alterations in mitochondrial activity play a significant role in the aging process and the 
progression of stroke, cancer, and neurodegenerative diseases such as Alzheimer's Disease 
[2—4]. Monitoring these alterations could be useful as a biomarker for characterizing disease 
progression or determining the efficacy of novel therapies. Demand is therefore high to 
develop robust tools for monitoring metabolic activity with high spatial and temporal 
resolution in vivo, in which the complex interactions between neurons, astrocytes, and 
cerebral vasculature remain preserved [5]. Reduced nicotinamide adenine dinucleotide 
(NADH) is the ubiquitous electron carrier that plays crucial roles in both glycolysis and 
oxidative metabolism, and monitoring its autofluorescence has proven useful as a minimally 
invasive indicator of metabolic activity for -50 years [6-8]. With absorption and emission 
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peaks in the UV region, the utility of monitoring NADH in scattering tissue samples with 
high resolution remained limited until the development of 2-Photon (2P) microscopy-based 
techniques, which capitalize on the spatially confined nonlinear excitation effect and the deep 
penetrating nature of near infrared excitation light [9-13]. More recent methods have coupled 
2P microscopy with time -resolved fluorescence lifetime imaging microscopy (FLIM), 
extending the utility and specificity of NADH autofluorescence measurement to better 
characterize microenvironmental conditions around the fluorophore [14-17]. Using time- 
correlated single photon counting (TCSPC) methods, novel FLIM and anisotropy studies of 
NADH model its fluorescence decay as the sum of two or more decaying exponentials, with 
each term reportedly representing a different enzymatic bound formulation, or "specie" 
[18,19], of NADH. These species ostensibly could be more specific biomarkers of molecular 
pathways associated with glycolysis or oxidative phosphorylation [18-23], but more 
extensive characterization is required. Here, we describe 2P- FLIM based in vivo observations 
of cerebral NADH in the cortices of anesthetized rats using a customized imaging system and 
analysis software. We observed four distinct cerebral NADH species under physiological 
conditions with high spatial resolution, and we monitored transient alterations in the relative 
amounts of these species under anoxic conditions with relatively high temporal resolution. 

2. Methods 

2. 1. Customized imaging system 

Figure la depicts the 2P FLIM portion of our multimodal microscope [24,25], custom- 
designed for in vivo investigations of cerebral metabolism in small animals. Excitation is 
provided by a tunable Ti:sapphire laser (Mai Tai, Spectra Physics), delivering -360 fs pulses 
after the objective at 80 MHz and tuned to X = 740 nm. The excitation intensity is modulated 
with an electro-optic modulator (ConOptics, Danbury CT), raster-scanned across a field of 
view (FOV) ranging between 75 and 300 (xm using galvanometer-based scanner mirrors 
(Cambridge Technology, Inc), and focused on the sample by a water immersion objective 
(Olympus XLumPlan Fluor, 20X, 1.00 NA, 2 mm working distance). Using a series of 
dichroic and bandpass filters, emitted light is spectrally separated and detected by a custom- 
designed four-channel detector array with high collection efficiency using photomultiplier 
tubes (PMTs) (Channel 1: 680 ± 30 nm, Channel 2: 595 ± 25 nm, Channel 3: 525 ± 25 nm, 
Channel 4: 460 ± 30 nm). NADH fluorescence is collected in channel 4 by a novel hybrid 
PMT with high detection efficiency and minimal afterpulsing (HPM-100-40 Becker & Hickl 
GmbH). As detailed in Subsection 2.4, fluorescence signal of the exogenous dye 
sulfurhodamine 101 (SR101) was detected in channel 2 simultaneously with the NADH 
collection. The SRI 01 intensity was used to co-register the NADH image with cellular 
structures and to correct for absorption changes associated with altered cerebral blood flow 
during physiological manipulations. 

For FLIM investigations described herein, scanning and data collection were performed 
using commercial time-correlated single photon counting (TCSPC) hardware (SPC-150, 
GVD-120, DCC-100, Becker & Hickl GmbH) and control software (SPCM, Becker & Hickl 
GmbH). 

NADH autofluorescence is spectrally identical to that of nicotinamide adenine 
dinucleotide phosphate (NADPH). In principle, the detected autofluorescence could arise 
from a mix of both NADH and NADPH. In metabolically active brain tissue, however, 
NADPH is widely believed to contribute minimally to the fluorescence signal, due to its 
much lower concentration in brain tissue [26,27], low quantum yield [28], and insensitivity to 
metabolic perturbation [13,19,29]. 
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Fig. 1 . (a) 2-photon imaging portion of our custom built in vivo imaging system, modified for 
FLIM measurements. SH: shutter, M: reflecting mirror, P: polarizer, EOM: electro-optic 
modulator, XY: galvanometer-based scanners, DM: Dichroic Mirror, PCM: photon counting 
module, PMT: photomultiplier tube, (a) Inset: reflectance image of the sealed cranial window 
model, approximately 3 mm in diameter at 570 nm. (b) In vivo image of NADH 
autofluorescence in a cortical astrocyte (Scale bar: 1 0 um), with (c) pixel- wide distribution of 
2-component fit performed with commercial SPCImage software, clearly showing 4 distinct 
peaks, (d) Example 4-component lifetime fit of in vivo cortical NADH fluorescence decay, 
computed with custom software. Blue profile: Collected, binned photon counts, Black-dashed 
lines: fitting boundaries, Green profile: computed IRF, Red profile: Computed decay profile 
(e) residuals for computed 2-, 3-, and 4-component lifetime fits to the example in vivo lifetime 
data, with reduced yj errors of 195.16, 2.64, and 2.54, respectively. 

2.2. Imaging protocol 

High resolution (256 x 256 pixels) images were first acquired under baseline physiological 
conditions using At «50 ps binning intervals and 256 temporal channels (12.8 ns-long decay 
profiles) over a field of view (FOV) ranging from 75 to 200 Ltm. Excitation intensity was 
adjusted to yield photon count rates of approximately 500,000 counts per second for these 
measurements, as indicated by the count rates continuously displayed in the SPCM software. 
Although incident power varied with cortical depth, it remained well below 50 mW for all 
measurements. For a given baseline measurement, such as the one provided in Fig. lb, a 
repetitive raster scan was performed over the FOV at frame intervals of up to -900 ms for 120 
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s. The FLIM system accumulates photons over all frames into time-resolved photon 
distributions, one distribution for each pixel, yielding a 256 x 256 x 256 data cube resolved in 
x, y, and t with high signal to noise ratio. We performed an additional 3x3 spatial pixel 
binning to obtain decay profiles with -5000 photons. We regard this value as the approximate 
minimum number of photons required to perform our multiple-component lifetime fits, 
containing 2 constrained lifetimes and 2 free lifetimes and described in Subsections 2.3 and 

3.2, at each pixel [30]. 

After collecting physiological baseline measurements, a time series of measurements over 
the same FOV was collected with 32 x 32 pixel resolution. Spatial resolution was sacrificed to 
decrease the time-point measurement interval for the time series to 2.5 s, enabling time- 
resolved FLIM measurements during metabolic manipulations. For each time -point 
measurement, the FOV was repeatedly raster scanned for 2.5 s. Each individual frame 
required -50 ms. The software accumulated photons over all -50 frames and temporally 
binned using the same At « 50 ps intervals. Each time series consisted of continuous time- 
point measurements acquired for 5 minutes. To minimize photobleaching, excitation intensity 
was reduced to yield photon count rates of approximately 50,000 counts per second for these 
time -resolved measurements. Fluorescence decay measurements were added over all pixels in 
the field of view. For a given time -point measurement, decay profiles were added with their 
two neighboring time points, yielding temporally smooth transient profiles. Each decay 
profile for each time point was comprised of 400,000 recorded photons. The total greatly 
exceeds the necessary amount for our fitting procedure, yielding more robust and temporally 
smooth profiles. These experiments constitute our global characterization of cerebral tissue's 
response to anoxia. Subsequent investigations will feature modified experimental protocols 
with pointwise measurements, enabling both high spatial and temporal resolution for 
characterizing more localized, cell-specific changes in cerebral metabolism. 

2.3. Data analysis: instrument response function and lifetime fitting 

As an initial step for each experiment, the lifetime calculation accuracy was verified, and 
simultaneously, the instrument response function (IRF) was computed, using FLIM 
measurements of dissolved NADH (0.5 mM in saline). Early attempts to characterize the IRF 
using second harmonic generation (SHG) of collagen and sucrose were confounded by 
mismatch in spectral properties [17,31]. The IRF in a two-photon microscope can be recorded 
by using SHG or Hyper-Rayleigh scattering [32]. Both methods deliver ultrafast signals at a 
wavelength detectable by the FLIM system. However, in contrast to fluorescence, most of the 
light is emitted in forward direction. The use of SHG may also require a change in laser 
wavelength. This can cause a shift in amplitude and phase of the timing reference pulses and 
thus a shift in the measured IRF. As NADH lifetime decay is well characterized [19,21] and 
its spectral features well match cytosolic and enzymatic bound NADH in vivo, NADH 
solution serves as a suitable reference emitter for evaluating these system characteristics 
under experimental conditions identical to our in vivo baseline measurements. 

The IRF was computed using an iterative procedure. No a priori assumptions were made 
for the NADH decay parameters. For the first iteration, the decay portion of the experimental 
NADH solution measurement (from the peak time to the final time point, \f lna i) was fit as the 
sum of two decaying exponentials and a constant offset. A preliminary IRF was then 
computed by the expression 



where FFT and FFT -1 denote the Fast Fourier Transform and inverse Fast Fourier Transform, 
respectively. For each subsequent iteration the entire experimental decay curve of measured 
NADH solution (from t = 0 to tf ma i) and the computed IRF{t) from the previous iteration were 




(1) 
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used to compute coefficients for a new theoretical decay. The new theoretical profile was then 
used in Eq. (1) to re -compute IRF{t). The procedure continued for 5 iterations. The decay 
coefficients yielded consistent values after the second iteration. The coefficients and IRF(t) 
from the iteration yielding the smallest fitting error (reduced x 2 ) were chosen. 

The fluorescence decays measured from in vivo brain tissue were fit with the expressions 



N 

I ,heore ti cal( t ) = I offse, +J E a i eX V 



(2) 



Kom (0 = imt) ® ( v (0 + 4— (0) • ( 3 ) 

where N = 4, as described in Subsection 3.2, a, and r, represent the amplitude and lifetime of 
component 7, and I 0 ff set represents the constant baseline intensity resultant from detector noise 
or ambient room light. I 0 A(t) models the ultrafast signal component observed in the detected 
response. This component generally constitutes a miniscule amount of the detected response 
and its origin is not well understood. In our experiments, the ultrafast component likely arises 
from backscattered excitation light [19]. Generally, negligible amounts of excitation light 
penetrate through the emission filters and reach the detectors. However, in highly scattering 
tissue such as the brain, the contribution of the ultrashort excitation pulse should be accounted 
for and can be modeled as a Dirac delta function [19,33]. 

2.4 In vivo animal studies 

Animals were prepared under a protocol approved by the Subcommittee on Research and 
Animal Care at Massachusetts General Hospital. As described previously [34] and illustrated 
in Fig. la (inset), cerebral cortices were exposed via a sealed cranial window in the skulls of 
anesthetized and ventilated Sprague Dawley rats (male, 195-265 g, Charles River). Cortical 
astrocytes were fluorescently labeled by topical application of sulforhodamine 101 dye 
(SR101, Sigma) [35]. Rats were initially anesthetized by isoflurane (1-2%). Following 
preparation under isoflurane, rats remained anesthetized and immobilized for the duration of 
the experiment with continuous infusion of alpha-chloralose (25 mg/kg/hr, Sigma Aldrich) 
and pancuronium bromide (2 mg/kg/hr, Hospira), respectively. Body temperature was 
monitored and controlled using a rectal temperature probe and heating blanket, respectively. 
The right femoral artery and vein were cannulated to monitor blood gases and pressure and to 
administer anesthetic. Using a flowmeter and a ventilator, gas fractions and inspiration rates 
were adjusted to maintain systemic blood gas and pH levels within physiologically normal 
range (pC0 2 : 37 ± 5 mmHg, p0 2 : 115 ± 15 mmHg, pH 7.39 ± 0.05). Brief periods of anoxia, 
not exceeding 45 s, were administered by respiratory arrest. 

Reduction of inspired oxygen stimulates increases in cerebral blood flow and volume in 
vivo, resulting in a greater amount of hemoglobin throughout the cerebral vasculature. We 
corrected for the adverse influence of augmented hemoglobin absorption and obstructed 
excitation light using the SRI 01 transient profile and a correction algorithm published 
previously [34]. Only the SR101 intensity signal was required for the correction algorithm. 
For each time point, the intensity was computed as time integrated SRI 01 photon counts at 
each pixel. No FLIM fitting was performed on the SR101 data. To verify that SR101 does not 
interfere with NADH lifetime measurements, we imaged solutions of SR101 in saline. 
Although the detected SRI 01 intensity in channel 2 was comparable to that observed in our in 
vivo experiments, negligible amounts of SR101 fluorescence were detected in the NADH 
channel from SRI 01 solutions. 
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3. Results 



3.1. FLIM measurements of NADH solution 

The fluorescence decay profile of a pure NADH solution was modeled as a bi-exponential 
decay, in which the two terms reportedly account for different folding conformations of the 
molecule [19,36]. Using a nonlinear least squares algorithm [15], we consistently computed 
fluorescence lifetimes of 395 ± 10 ps and 1.05 ± 0.010 ns for the two components of NADH 
solution, with relative amplitudes of 90% and 10%, respectively. Though our observations 
yield an average lifetime similar to previous 2P microscopy characterizations of NADH 
solution (460 ps, compared to 444 ps [19] and 480 ps [21]), the computed lifetimes of the 
short and long components differ by 12.5% and 24% of reported results. The discrepancy may 
be attributable to different micro environmental conditions such as solvent medium, 
temperature, pH, viscosity, or concentration. 

3.2. In vivo measurements: baseline conditions 

Figure lb displays an autofluorescence intensity image of a cortical astrocyte in vivo i 
measured approximately 65 (xm below the tissue surface. The astrocytes' processes appear 
relatively bright; however, the image appears diffuse and noisy due to the nonzero signal 
outside the cell. The extracellular intensity arises from NADH autofluorescence in the 
surrounding neuronal and astrocytic processes. Analyses of preliminary in vivo brain tissue 
experiments were performed with commercial software (SPCImage, Becker & Hickl GmbH). 
Attempts to fit baseline experimental in vivo data at each pixel (after 3x3 binning) as the sum 
of two or three decaying exponentials routinely indicated the presence of four distinct NADH 
components, illustrated in Fig. lc and consistent with previous reports [19,37] that detail in 
vitro brain slice measurements and rat liver, as well as our own unpublished observations in 
brain slices. Each of the four components were found to be independently sensitive to 
metabolic perturbations such as varied fraction of inspired oxygen or application of potassium 
cyanide [unpublished results]. This motivated the development of custom lifetime-fitting 
software in Matlab for all subsequent data analysis. 

From our early fitting computations, which computed all lifetimes and amplitudes, we 
observed that the first two components consistently yielded lifetime values close to 0.395 and 
1.0 ns. This lead to the conclusion that the first two components (a) correspond to the 
different folding conformations of free NADH identified from our NADH solution 
measurements and (b) the lifetimes of these components remained the same between our 
NADH measurements in solution and those measured in vivo. Consequently, we modified our 
fitting routine such that Tj and r 2 were constrained within ±10 ps of the lifetime values 
computed for NADH solution at the start of each experiment. A nonlinear least squares fitting 
routine was then used to compute I 0 , r 3 , r 4 , and a, for i = 1-4. Figure Id illustrates example 
profiles of in vivo cerebral NADH fluorescence, including the I measure d (t), hwdei (t), and IRF{t), 
and the associated fit values computed with our custom software. Residual values multi- 
exponential fits containing 2, 3, and 4 components are displayed in Fig. le. A bi-exponential 
model was found to be inadequate to model in vivo NADH fluorescence lifetime, as indicated 
by the residual profile. While both 3- and 4-component models yielded more satisfactory 
residual profiles, the 4-component model consistently resulted in smaller fitting errors. 

Figures 2a-2b display representative high-resolution intensity images of rat cortical tissue 
collected in vivo with our microscope, roughly 70 (xm below the brain surface. Astrocytic cell 
bodies appear bright in the SR101 fluorescence images due to their preferential dye uptake 
over neurons [35]. Their processes envelop cerebral microvasculature, whose cross sections 
appear as dark circles in the intensity images. As shown in Fig. 2a, the SRI 01 intensity image 
permitted identification of astrocytic cell bodies (yellow) and blood vessels (red). We 
designated the remaining tissue as the neuropil, comprised primarily of processes from both 
neuronal dendrites and astrocytes. While select cell bodies appear brighter than the 
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Fig. 2. Spatial distributions of cerebral NADH components in vivo (a-b) Example in vivo 
intensity images collected from the rat cortex, ~70 urn below the cortical surface, (a) SRI 01 
fluorescence preferentially accumulates within astrocytic cell bodies and processes, enabling 
distinction between astrocytic cell bodies (yellow ROIs), blood vessels (red ROIs), and 
neuropil, (b) NADH fluorescence can be seen more uniformly throughout the tissue, though 
overlying vasculature results in shadows, (c) Distributions of fractional fluorescence for the 
four computed NADH components and purported scattering component, determined by fitting 
the fluorescence decays at each pixel and classifying pixels as astrocytes (yellow ROIs) blood 
vessels (red ROIs), or surrounding neuropil (scale bar: 50 um, error bars: standard error) 



background, NADH autofluorescence is generally ubiquitous within the tissue, allowing for 
little contrast of structural features at this magnification. As seen in the lower left corner, 
shadows of overlying blood vessels are visible as a result of obstructed focusing of excitation 
light by scattering red blood cells and from absorption of NADH emission by hemoglobin 
[34]. Figure 2c displays the tissue-specific fractional fluorescence of each NADH component, 
calculated at each pixel as 



Figure 2c also shows the average amount of detected ultrafast signal in each tissue (7 0 in 
Eq. (3), which we attribute to back-scattered excitation light, normalized to its overall 
maximal value. We observed that the NADH components generally appear heterogeneously 
distributed and diffuse across the tissue, and the components' spatial distributions do not 
coregister well with particular cellular structures seen in the SRI 01 intensity images. This 
phenomenon was observed in fractional fluorescence images computed with both our custom 
software and with commercial SPCImage software, used to calculate 2- and 3 -component fits 




(4) 
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[unpublished results]. The first NADH component (CI), which has the shortest lifetime (fixed 
at -0.4 ns) and likely represents NADH unbound from any enzyme, constitutes a similar 
amount of the fractional fluorescence for astrocytes, blood vessels, and the neuropil. Its 
relative amplitude suggests that it is the most prevalent component. However, CI makes the 
smallest contribution to the overall fluorescence signal due to its shorter lifetime. The lifetime 
of C2 was constrained to ~1 ns. Its computed relative amplitude was found to be less than that 
of CI but greater than the other components. Its contribution to the fractional fluorescence 
was smaller in blood vessels compared to astrocytes and neuropil. C3 has a computed peak 
lifetime of 1.7 ns. This lifetime value corresponds well with a component observed in 
solutions of NADH bound to both malate dehydrogenase and lactate dehydrogenase [20,38]. 
C3 appears in similar amounts as C2 in both astrocytes and the neuropil, but it contributes 
considerably more than C2 to fractional fluorescence in the blood vessels. C4 has a peak 
lifetime value of 3.2 ns. In astrocytes, blood vessels, and neuropil, it has the greatest 
contribution to the fractional fluorescence, likely attributable to its longer lifetime. A 
component with similar lifetime and relative proportion has been observed in malignant and 
normal HS578 breast cells [20]. The contribution from backscattered excitation light is 
substantially higher in astrocytic cell bodies than the neuropil, and it is higher still in the 
blood vessels. This suggests that the back scattering in tissue arises primarily from cell bodies 
rather than dendritic processes. 

Our baseline measurements consisted of multi-exponential lifetime fitting at each pixel 
using fluorescence profiles consisting of approximately 5000 photons. We consider 5000 to 
be the minimum number of photons required to perform rough lifetime fits using our method 
of constraining the first 2 of 4 lifetime values. Table 1 displays our modeling results for a 
theoretical decay consisting of 5000 photons with various added amounts of noise. The 
amplitude of the noise ranged between 0 and 200% of shot noise (the square root of the 
photon counts in each binning interval), and ii and x 2 were constrained within ± 10 ps. The 
mean and standard deviations are provided for the coefficients from a total of 10 trials for 
each noise amplitude. 

In the presence of shot noise, the computed lifetime values were found to be reasonably 
stable, while computed amplitudes demonstrated more variability. Average lifetime values 
remained within 10% of the original value, and amplitudes varying by as much as -35%. 
Adding twice as much noise yielded more variability, with up to 45% change in computed 
lifetime value and 65% change in amplitude. Under normal conditions, the noise component 
should remain well below twice the shot noise. Therefore, we believe analysis of profiles with 
5000 photons is suitable for measuring rough but reasonably stable FLIM parameters. 



Table 1. Computed parameters from profiles modeled with varied levels of added noise 





Original 


+ 0.5x shot noise 


+ lx shot noise 


+ 2x shot noise 




73.42 


72.57 ±5.79 


69.85 ± 6.83 


76.08 ± 9.27 


a 2 


64.79 


58.54 ±20.81 


69.79 ± 9.04 


61.84 ±26.50 


<*3 


18.22 


26.23 ± 17.01 


20.99 ±6.86 


35.41 ±22.40 


a 4 


27.07 


29.32 ±3.09 


28.82 ±2.59 


23.66 ± 11.58 


Tl 


0.41 


0.399 ±0.01 


0.40 ± 0.01 


0.41 ±0.01 


Tl 


0.99 


1.00 ±0.01 


0.99 ±0.00 


0.99 ±0.01 


T} 


1.82 


1.51 ±0.23 


1.67 ±0.22 


2.21 ±0.90 


T 4 


4.21 


4.37 ±0.22 


4.62 ±0.15 


6.02 ±2.12 


Itiffset 


0.05 


0.04 ±0.01 


0.04 ± 0.02 


0.04 ± 0.02 



3.3. In vivo measurements: anoxia transients 



Figure 3 a displays representative transient profiles of amplitude weighted lifetimes, a,r„ i: 1- 
4,; of the four NADH components during a brief period of anoxia, induced by respiratory 
arrest. The profiles were computed from decay profiles assembled over the entire field of 
view. The profiles have also been normalized to the average baseline value and corrected for 
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Fig. 3. NADH transients during anoxia (a) Representative global transients of NADH 
components (amplitude-weighted lifetime, a,r,), and intensity transients for NADH, SR101, 
and backscattering. The shaded region corresponds to a 45 s period of anoxia, after which 
respiration was immediately restored, (b-e) Average component-specific changes in (b) 
anoxia-induced onset time (c) rise time (d) relative undershoot, and (e) maximal relative 
change, as measured from 10 trials of anoxia. Error bars denote standard error. * denote 
statistically significant differences from indicated components, measured with paired t-tests (a 
= 0.05) 

changes in blood volume and hemoglobin using the metabolically inert SRI 01 dye signal and 
a correction algorithm [34]. NADH intensity, SR101 intensity, and Iq are also included. 
During the onset of anoxia, pronounced, steep rises are observed in all NADH components, 
NADH intensity, and /o(0 as halted oxygen supply progressively limits oxidative 
phosphorylation. Transient profiles were analyzed for 10 anoxia trials in a total of 7 rats, and 
results are provided in Figs. 3b-3e. In all trials, respiration was restored after 45 seconds. 
Arterial blood pressure dropped drastically from -110 mmHg to ~50 mmHg during the 
anoxic episode, but quickly recovered during oxygen restoration. Onset time (f 10 ) and rise 



#179359 - $15.00 USD Received 6 Nov 2012; revised 17 Jan 2013; accepted 17 Jan 2013; published 22 Jan 2013 
(C) 2013 OSA 1 February 2013 / Vol. 4, No. 2 / BIOMEDICAL OPTICS EXPRESS 317 



time (tgo-tio) were calculated from interpolated time courses (At = 0.1 s) by rinding the times 
t w and t go required to increase from average baseline value to 10% and 90% of the maximum 
value, respectively. For all measured parameters, onset times occurred approximately 15 s 
after initiating respiratory arrest. Paired ?-tests (a = 0.05) were performed for all 
measurements in all animals. Onset times and peak undershoot values were found not to be 
significantly different across the 4 NADH components or I 0 . NADH components 1 and 3 had 
significantly different rise times, with component 3 rising nearly 3 times faster than 
component 1. Components 1, 2, and 3 all have significantly different maximal changes, with 
component 1 changing by nearly 2 and 3 times as much as components 2 and 3, respectively. 
I 0 represents the changes in magnitude of backscattered excitation light. Under hypoxic 
conditions, reports have suggested that increases in scattering can be attributed to 
mitochondrial swelling [19,39]. The observed transient profiles in the backscatter signal were 
found to correlate well with profiles of NADH component 4, with similar onset times, 
amplitude changes, and undershoots. While more rigorous investigation is required for 
validation, the correlation suggests that NADH component 4 may represent a formulation of 
NADH bound to a large mitochondrial enzyme, and could be a more specific indicator of 
changes in oxidative metabolism. 

4. Discussion and conclusion 

NADH fluorescence lifetime has been characterized in a wide range of biological samples 
and pathological conditions, including isolated mitochondria [33,37], in vitro cell cultures 
[20,22,29,37,39-41], ex vivo tissue [19,21,23,37,42], and in vivo measurements [23,43,44]. 
To our knowledge, we report the first observations of NADH lifetime using 2P microscopy in 
cerebral tissue in vivo. NADH fluorescence decays are typically modeled as the sum of 
multiple (2 - 4) decaying exponentials, representing different enzyme -bound formulations, 
(i.e. species [19]), with distinct lifetimes. Although the species with shortest lifetime is widely 
believed to indicate 'free' NADH, unbound from any enzyme, the total number of identifiable 
species, and their associated lifetimes and represented enzymatic processes, remains 
unresolved. NADH participates in numerous reactions of glycolysis, the Kreb's cycle, and 
oxidative phosphorylation. Several enzymes and proteins exist to which NADH can 
intermediately bind both in the cytosol and mitochondria. Therefore, practical models of 
multiple exponential decays necessarily represent an oversimplification of the underlying 
cellular metabolic processes [45]. Investigators have reported similarities in lifetime values 
between their intracellular NADH measurements and solutions of NADH complexed with 
isolated enzymes such as lactate dehydrogenase and malate dehydrogenase [19,20]. Together 
with mitochondrial complex I, it is likely that the observed NADH species reflect 
intermediate binding with these particular enzymes. However, more investigation is required, 
particularly to account for the differences in enzymatic activity and environmental conditions 
between solution medium and the intracellular environment. 

Discrepancies in the number of components and lifetime values exist between our 
observations and published values. As one of several examples, Niesner et al. reported a 
distribution of NADH lifetimes suggesting the presence of only 2 components, representing 
unbound and "protein-bound" NADH, observed in cultured dermal fibroblasts [18]. The 
discrepancies between our observations and those of Niesner et al. and others may be 
attributable to several different experimental parameters and conditions utilized in each study. 
In particular, the different cell types, investigated primarily under in vitro and ex vivo 
conditions, and different computational methods would likely account for most of the 
discrepancies. To our knowledge, we report the first in vivo FLIM observations of NADH in 
the rodent cerebral cortex. Under these in vivo conditions, we monitored the metabolic 
activity of more than one cell type (neurons, astrocytes, vascular endothelial cells, etc), and 
the supply of oxygen and other metabolites are likely more spatially heterogenous than under 
cell culture or brain slice conditions [24,46]. Niesner et al. observed that the weighted 
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averages of NADH lifetimes in dermal fibroblasts differ from those measured in solution, and 
they suggested that the differences are attributable to the variations in microenvironmental 
conditions such as viscosity, pH, and refractive index. These environmental differences would 
likely be more pronounced and heterogeneous under in vivo conditions, and this may explain 
why our lifetime distribution in Fig. lc differs markedly from their distribution. Additionally, 
they utilized the noniterative Prony method rather than the nonlinear least squares method to 
calculate their lifetimes, which also could contribute to the difference. Our observation of 4 
NADH components is consistent with the results of Vishwasrao et al. [19], but it differs from 
those of Chia et al. [21]. Both investigators measured NADH lifetime in hippocampal slices 
of Sprague Dawley rats. Compared to our results, Chia et al. performed analysis on traces 
with fewer photon counts (-1200), and utilized the commercial SPCIMage software, which 
can perform fits with a maximum of 3 components. This may account for the differences in 
our results. 

Like NADH, another important electron carrier for oxidative metabolism, flavin adenine 
dinucleotide (FAD), is also autofluorescent, but in its oxidized form. In principle increased 
FAD fluorescence suggests an increase in oxidative phosphorylation via succinate 
dehydrogenase (mitochondrial complex II) [7]. Simultaneous detection of reduced NADH 
and oxidized FAD fluorescence enables calculation of the redox ratio, a semi-quantitative 
measure of energy respiration and the oxidation-reduction state in the mitochondrial matrix 
that is independent of the factors influencing the fluorescence measurements [22,47—49]. 
FAD measurements are reportedly limited by low signal, spectral overlap with other 
endogenous fluorophores such as lipofuscin [50], and less responsiveness to metabolic 
perturbations [51]. Efforts are currently underway to address these limitations and collect 
robust 2P in vivo FLIM measurements of FAD in channel 2 of our system. 

Our in vivo observations suggest the existence of a minimum of 4 distinct, detectable 
species in brain tissue. This motivated us to evaluate the lifetime data using a 'quasi-global' 
analysis method, where profiles were fit to a 4-component multi-exponential fit, with 2 
lifetimes held constrained to those measured in NADH solution, our designated reference 
emitter [17,52]. Further, rigorous investigation is required to identify which metabolic 
enzymes and proteins could be represented by each specie. In order to monitor shifts in 
glycolytic and oxidative metabolism, however, this level of detail may not be necessary. It 
may be sufficient to associate the species with different metabolic reactions without 
identifying its associated enzyme. Future investigations will explore the viability of this 
technique by monitoring transient features of each component during exposure to various 
metabolic inhibitors. 

Our identification of 4 distinct components agrees well with reported ex vivo brain tissue 
and isolated liver mitochondria measurements [19,37]. Though the physical significance of 
parameters derived from a 4-component exponential fit has been disputed [53], our 
hypothesis is supported by the observation that each component responds differently to 
metabolic perturbations. For our baseline in vivo measurements, we performed rough multi- 
exponential lifetime fits with high spatial resolution. Fits were performed at each pixel on 
profiles with -5000 photons, after 3x3 binning. Modeling results also showed that calculated 
lifetime parameters were reasonably stable for simulated fluorescence decays with -5000 
photons and added shot noise. Published values for the minimum number of photons required 
for precise resolution of 2 free lifetime components range from 1000 to 400,000 [30,54]. 
Although Kollner and Wolfram theoretically calculated that 400,000 photons are needed to 
determine the components of a double-exponential decay [54], it should be noted that they 
used an example of an extremely unfavorable composition of the decay (10% of 2 ns and 90% 
of 4 ns). In practice, the chance of resolving a multi-exponential decay function into its 
components dramatically increases with the ratios of the component lifetimes, and with the 
amplitude(s) of the fast components. For the typical decay composition of NADH, we 
obtained reasonably stable parameters by fitting decay functions containing approximately 
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5000 photons. We therefore consider 5000 as the minimum number of photons required to 
perform our multiple-component lifetime fits, containing 2 constrained lifetimes and 2 free 
lifetimes. This minimum photon threshold, coupled with 3x3 spatial binning, was 
determined to be the best compromise between spatial resolution and SNR, and yielded 
consistent lifetime values across measurements. Nevertheless, the low 2P cross section of 
NADH [12] limits signal for 2P FLIM measurements using TCSPC. This, in turn, reduces the 
precision for performing multi -exponential fits for each pixel and limits the ability to 
precisely measure the distributions of multiple NADH components in cytosol and individual 
mitochondria in vivo. In principle, collecting more photons, using either higher excitation 
intensity, longer dwell times, or longer scanning duration could enable fits with greater 
certainty [53]. However, great care must be taken to prevent photobleaching. 

In contrast to previous reports [20,21], we did not observe punctate clusters of 
fluorescence in our lifetime images. Though their spatial distributions were heterogeneously 
distributed, fluorescence signals of each component could generally not be registered with 
distinct cellular structures. Instead, the fractional fluorescence images for our calculated 
components appear diffusely distributed. Our results likely differ from previous reports 
because we imaged with lower spatial resolution and a larger FOV. Imaging with a larger 
FOV is necessary for evaluating the complex metabolic interactions between different cell 
types associated with neurovascular coupling. Additionally, punctate intracellular features 
could also be blurred by potential motion artifacts associated with respiration during in vivo 
imaging. 

The amount of fractional fluorescence intensity from each NADH component differs 
subtly between astrocytic cell bodies, blood vessels, and neuropil, as does the amount of 
back-scattered excitation light. C3 and C4 are believed to represent NADH bound to 
mitochondrial enzymes, and they contribute most strongly to fluorescence in all tissue types. 
Erythrocytes are the primary cell type in blood and lack mitochondria. The prevalence of C3 
and C4 in blood vessels is therefore an unexpected observation. This finding requires more 
investigation, but the detected C3 and C4 could have been measured from other more 
metabolically active blood cells that contain mitochondria, such as macrophages and platelets. 
Blood vessels also scatter considerably more excitation light than astrocytes or neuropil. This 
could imply that the erythrocyte membranes have higher scattering properties than brain 
tissue. 

Due to its spectral compatibility, measuring a solution of NADH and fitting to a 2- 
component model proved to be a more reliable method to determine the IRF compared to 
SHG measurement of collagen or sucrose. We employed methods of the global analysis 
technique [16] by constraining the first 2 lifetime values in our in vivo data fits, yielding 
robust 4-component fits of NADH lifetime in the cerebral cortex. We believe that 
constraining the values of the first two components is intuitive and practical. While other 
microenvironmental factors such as temperature and viscosity influence the lifetime value, it 
is reasonable to expect that enzymatic binding has the biggest influence on NADH. 
Therefore, in our fitting routine, restricting the upper and lower bounds of components 1 and 
2 within 10 ps of the measurement grants the flexibility for these other factors while asserting 
that they correspond to distinct conformations of free NADH. 

For our measurements of in vivo metabolic transients, we performed precise multi- 
exponential lifetime fits with lower spatial resolution, a suitable approach for investigating 
global phenomena such as anoxia. Under normal physiological conditions, the brain is 
supplied with a surplus of oxygen. Our observations show that, under our experimental 
conditions, only —15 s of respiratory arrest is sufficient to deplete vascular oxygen supply to 
levels where oxidative phosphorylation rates start to diminish in the rat cerebral cortex. After 
this point, all observable NADH species increase sharply, as does tissue backscattering. The 
increase in scattering is believed to arise from mitochondrial swelling. However, it could also 
be attributable in part to anoxia-induced increases blood flow and hematocrit. Onset times for 
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these changes were found to be similar for all species and backscattering. However, the 
shortest lifetime component, which is believed to indicate unbound NADH, has the most 
pronounced increase. Kinetic features of the longest lifetime component are strikingly similar 
to those of tissue scattering, suggesting that this enzymatic formulation of NADH may be 
involved in the same metabolic process that induces mitochondrial swelling. These analyses 
over the entire FOV motivate further modifications to our technique to enable spatially 
resolved measurements with high temporal resolution. Future investigations will include 
point-wise measurements of fluorescence lifetime and anisotropy at designated locations in 
the field of view, enabling distinction of cell-specific responses to metabolic perturbations in 
the cortex. 
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